Dirac Equation and Quantum Relativistic Effects in a Single Trapped Ion 
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We present a method of simulating the Dirac equation in 3+1 dimensions for a free spin- 1/2 
particle in a single trapped ion. The Dirac bispinor is represented by four ionic internal states, and 
position and momentum of the Dirac particle are associated with the respective ionic variables. We 
show also how to simulate the simplified 1+1 case, requiring the manipulation of only two internal 
levels and one motional degree of freedom. Moreover, we study relevant quantum-relativistic effects, 
like the Zitterbewegung and Klein's paradox, the transition from massless to massive fermions, and 
the relativistic and nonrelativistic limits, via the tuning of controllable experimental parameters. 
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The search for a fully relativistic Schrodinger equation 
gave rise to the Klein- Gordon and Dirac equations. P. A. 
M. Dirac looked for a Lorentz-covariant wave equation 
that is linear in spatial and time derivatives, expecting 
that the interpretation of the square wave function as a 
probability density holds. As a result, he obtained a fully 
covariant wave equation for a spin-1/2 massive particle, 
which incorporated ab initio the spin degree of freedom. 
It is known [l[ that the Dirac formalism describes ac- 
curately the spectrum of the hydrogen atom and that 
it plays a central role in quantum field theory, where 
creation and annihilation of particles are allowed. How- 
ever, the one-particle solutions of the Dirac equation in 
relativistic quantum mechanics predict some astonishing 
effects, like the Zitterbewegung and the Klein's paradox. 

In recent years, a growing interest has appeared re- 
garding simulations of relativistic effects in controllable 
physical systems. Some examples are the simulation of 
Unruh effect in trapped ions [2j], the Zitterbewegung for 
massive fermions in solid state physics !S| , and black- hole 
properties in the realm of Bose-Einstein condensates [3| . 
Moreover, the low-energy excitations of a nonrelativis- 
tic two-dimensional electron system in a single layer of 
graphite (graphene) are known to follow the Dirac- Weyl 
equations for massless relativistic particles [H, On 
the other hand, the fresh dialog between quantum infor- 
mation and special relativity has raised important issues 
concerning the quantum information content of Dirac 
bispinors under Lorentz transformations [7[. 

In this Letter, we propose the simulation of the Dirac 
equation for a free spin-1/2 particle in a single trapped 
ion. We show how to implement realistic interactions on 
four ionic internal levels, coupled to the motional degrees 
of freedom, so as to reproduce this fundamental quantum 
relativistic wave equation. We propose also the simula- 
tion of the Dirac equation in 1+1 dimensions, requiring 
only the control of two internal levels and one motional 
degree of freedom. We study some quantum-relativistic 



effects, like the Zitterbewegung and the Klein's paradox, 
in terms of measurable observables. Moreover, we dis- 
cuss the transition from massless to massive fermions, 
and from the relativistic to the nonrelativistic limit. Fi- 
nally, we describe a possible experimental scenario. 

We consider a single ion of mass M inside a Paul trap 
with frequencies v x , v y , and v z , where four metastable 
ionic internal states, |a),|6),|c), and |d), may be cou- 
pled pairwise to the center-of-mass (CM) motion in di- 
rections x, y, and z. We will make use of three stan- 
dard interactions in trapped-ion technology, allowing 
for the coherent control of the vibronic dynamics 
First, a carrier interaction consisting of a coherent driv- 
ing field acting resonantly on a pair of internal levels, 
while leaving untouched the motional degrees of free- 
dom. It can be described effectively by the Hamilto- 
nian H a = Ml(tj + e 1 ^ + <7~e~^), where a + and a~ are 
the raising and lowering ionic spin-1/2 operators, respec- 
tively, and f2 is the associated coupling strength. The 
phases and frequencies of the laser field could be ad- 
justed so as to produce H ax = hVl x o~ x , H a = h£l y a y , 
and H a _, = Ml z a z , where a x , o~ y , and o~ z , are atomic 
Pauli operators in the conventional directions x, y, and 
z. Second, a Jaynes-Cummings (JC) interaction, usu- 
ally called red-sideband excitation, consisting of a laser 
field acting resonantly on two internal levels and one 
of the vibrational CM modes. Typically, a resonant 
JC coupling induces an excitation in the internal lev- 
els while producing a deexcitation of the motional har- 
monic oscillator, and viceversa. The resonant JC Hamil- 
tonian can be written as H r = hrj(l(a + ae z ^ r +er~ e~ l ^ r ), 
where a and are the annihilation and creation op- 
erators associated with a motional degree of freedom. 
r\ = kyJh/2M v is the Lamb-Dicke parameter where 
k is the wave number of the driving field. Third, an anti- 
JC (AJC) interaction, consisting of a JC-like coupling 
tuned to the blue motional sideband with Hamiltonian 
H h = hr]n(cr + a 1 e l<t > h + <j~ae^ b ). In this case, an in- 
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ternal level excitation accompanies an excitation in the 
considered motional degree of freedom, and viceversa. 

All these interactions could be applied simultaneously 
and addressed to different pairs of internal levels cou- 
pled to different CM modes. For example, it is pos- 
sible to adjust field phases to implement a simultane- 
ous blue and red sideband excitation scheme to form the 
Hamiltonian — ihr] x (l x a x (a J. — a x ) — 2r\ x A, x Sl x a x p Xl 
with i{a\ — a x )/2 — A x p x /h. Here, A x := •Jh/2Mv x 
is the spread in position along the x-axis of the zero- 
point wavefunction and p x the corresponding dimen- 
sioned momentum operator. The physics of H?* can- 
not be described anymore by Rabi oscillations. In turn, 
it yields a conditional displacement in the motion de- 
pending on the internal state, producing the so-called 
Schrodinger cat states 

(sHH- By further manipulation of 
laser field directions and phases, we can also implement 
H^l = 2r] y A y n y a y Py and H%* = 2rj z A z tt z a x p z . This 
kind of interactions has already been produced in the 
lab, under resonant [n[ and dispersive conditions 

We define the wave vector associated with the four 
ionic internal levels as 



|*) := *» + * b \b) + * c 



*d\d) = 



(i) 



We may apply simultaneously different laser pulses, with 
proper directions and phases, rj = r/ x = r/ y — r/ z , A = 

^ r£ ~ " ^ qj ~ " ^ g - r£ <y g , r£ y 

f2 z , to compose the following Hamiltonian acting on 



H U = 2 V An(af + a b x c ) Px + 2r,An(af - a b y c )p y 



bd > 



+2nAn{a a x c - a° d ) Pz + Ml(a^ c + a 
We rewrite Eq. @ in the suitable matrix form 





2nAD,(a -p)+iMl 



2r]An(a ■ p)-ihft 




(2) 



(3) 



where each entry represents a 2 x 2 matrix. The asso- 
ciated Schrodinger equation, H-a\^f) = iHd\^>)/dt, per- 
forms the same dynamics as the Dirac equation in 3+1 
dimensions for a free spin-1/2 particle, where |$) repre- 
sents the four-component Dirac bispinor. This is easily 
seen if we express the Dirac equation 



'at 



ih-^r = H-Dip = ( C( 5 ' P + (3mc 2 )ip (4) 



in its "supersymmetric" representation [l[ 



c(a • p) — imc 2 

c(a ■ p) + imc 2 



(5) 



Here, the 4x4 matrix a :— (a x , a y , a z ) = off-diag(<7, a) 
is the velocity operator, [3 := off-diag(— ila, H2), and 



c := 2r]Afl , mc 



Ml, 



(6) 



are the speed of light and the electron rest energy, re- 
spectively. The notorious analogy between Eqs. ([3]) and 
([5]) shows that the quantum relativistic evolution of a 
spin-1/2 particle can be fully reproduced in a tabletop 
ion-trap experiment, allowing the study of otherwise un- 
accesible physical regimes and effects, as shown below. 

In the Dirac formalism, the spin-1/2 degree of freedom 
is incorporated ab initio. Moreover, the Dirac bispinor in 
Eq. (fT]) is built by components associated with positive 
and negative energies, £7 D = ±^p 2 c 2 + m 2 c 4 . This de- 
scription is the source of diverse controversial predictions, 
as the Zitterbewegung and the Klein's paradox. 

The Zitterbewegung is a known quantum-relativistic ef- 
fect consisting of a helicoidal motion of a free Dirac par- 
ticle, a natural consequence of the non-commutativity of 
its velocity operator components, cat, with i = x,y, z. It 
can be proved straightforwardly [l[ that the time evolu- 
tion of the position operator r — (x,y, z) in the Heisen- 
berg picture, following df/dt — [r, H-£>]/ih, reads 



f(t) = r(0) + 



t 



4n 2 A 2 il 2 p 
2nAQp\ ihrjAfl 



^ e 2iH Dt /h_ 1 y (7) 



Here, the first two terms on the r.h.s. account for the 
classical kinematics of a free particle, while the last os- 
cillating term is responsible for a transversal "quivering" 
motion. If we consider a bispinor state with a peaked mo- 
mentum around po, \*f?o) = |a)®exp[— (p— po) 2 /2a 2 ], the 
Zitterbewegung frequency associated with the measurable 
quantity (r(t)} can be estimated as 



wzb « 2\E D \/n = 2JAn 2 A 2 tt 2 P 2 /h 2 + n 2 



l^Nr] 2 n 2 



Q 2 



(8) 

where £b = (Hb) and N = (a^a) are the average en- 
ergy and phonon number, respectively. Similarly, we can 
estimate from Eq. (J7J) the Zitterbewegung amplitude as- 
sociated with (r(t)) as 



Rzb — 



2mc \ E 



r)h 2 VlVLA 



Ar] 2 n 2 A 2 p 2 + h 2 n 2 



(9) 



and Rzb ~ A, if nil ~ fi. 

The standard explanation of this erratic motion for a 
free Dirac particle invokes the interference between the 
positive- and negative-energy components of the Dirac 
bispinor following the dynamics in Eqs. © and ©. The 
predicted values for a real electron, luzb ~ 10 21 Hz and 
Rzb ~ 10 _3 A, are out of experimental reach, the effect 
has never been observed, and its existence is even ques- 
tioned by quantum field theory considerations. To sim- 
ulate quantum-relativistic effects in other physical sys- 
tems, like trapped ions or graphene, is not aimed at prov- 
ing their existence, but at exploiting the differences and 
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analogies in each field. Given the flexibility of trapped- 
ion systems, we will have access to a wide range of tun- 
able experimental parameters allowing for realistic 
and measurable uj^b ~ — 10 6 Hz and Rzb ~ — 10 3 A, 
depending on the initial vibronic states. Due to our piece- 
wise build-up of the 3+1 Dirac Hamiltonian of Eq. @, we 
can strongly reduce the experimental demands to study 
Dirac equations in 1+1 and 2+1 dimensions [3, [Hj]. In 
those cases, the Clifford algebra that characterizes the 
Dirac matrices is satisfied by the anticommuting Pauli 
matrices, {ai,<Tj} — 2Sij, where {A, B} := AB + BA, 
given that we only need 2(3) anticommuting matrices in 
the 1+1(2+1) case. Accordingly, the four components of 
the Dirac bispinor are conveniently reduced to only two. 

We focus now on the 1+1 dimensional case, which 
could be reached by current experiments, keeping most 
striking results available. After a unitary transforma- 
tion around the rc-axis, transforming a v into a z , the 1+1 
Dirac equation stemming from Eq. |2]) can be cast into 
ihd\¥^)/dt = JT^I^W) with 



H 



(i) 



ihrjO,cr x (a x — a x ) + hfla z 
2rjAfla x p x + hQa Zl 



(10) 



where the Dirac "spinor" = *^ 1J |a) + *^|6>. Note 

that these two components are not associated with the 
spin-1/2 degree of freedom but are linear combinations of 
positive- and negative-energy solutions Our setup 

is now reduced to a simultaneous JC + AJC interaction, 
a x p Xl and a Stark-shift term a z acting on two internal 
levels. In Eq. (JTDJ), the use of a z is for the sake of peda- 
gogy, and returning to a y will not affect the results. 

In the nonrelativistic limit, mc 2 3> p x c, or equivalently 
n 3> r/fl, we may identify the dispersive limit in Eq. (|10|) . 
and derive the squeezing second-order Hamiltonian 



(1) 2fAHl 2 pi 



(ii) 



yielding the expected Schrodingcr Hamiltonian associ- 
ated with the classical kinetic energy of a free particle. In 
the relativistic limit, mc 2 <C p x c, which includes m = 0, 
the 1+1 case reduces to HjP = ihr)Qcr x (a x — a^.which 
produces Schrodinger cats, as commented above [9L fTol| . 

For a massless particle, we can show that dx/dt = 
[x,H^ ] ]/ih = 2r)Atta x and da x /dt = [a x ,H^ ] ]/ih = 0. 
In consequence, a x is a constant of motion and the time 
evolution of the position operator, 



x(t) = x(0) + 2nAn<j x t , 



(12) 



is classical and does not involve Zitterbewegung. On 
the other hand, for a massive particle, da x /dt = 
[a x , Hi)}/ih = — £la y , doyjdt ^ 0, and so forth. The 
produced set of differential equations yields a Zitterbewe- 
gung oscillatory solution similar to Eq. ([7]) , with the evi- 
dent simplification via the replacements of by 



P by p x , and a by a x . It is noteworthy to mention 
that the phenomenon of mass acquisition, which could 
be done here in a continuous manner by raising the cou- 
pling strength Q, is related to the spontaneous symmetry 
breaking mechanism of the Higgs field. 

At this stage, it is clear that the measurement of the 
expectation value of the position operator, (x(t)), as a 
function of the interaction time t is of importance. A 
recent proposal for realizing fast measurements of mo- 
tional quadratures 14| relies on the possibility of measur- 
ing the population of an ionic internal level, P a , at short 
probe-motion interaction times r with high precision [l5| . 
Given an initial vibronic state p(0) = |+^)(+^| <8 p m , 
where |±^) = (\a) ± e % ^\b)) / 'y/2 and p m describes an un- 
known motional density operator, we can make use of 



'(r) 



(13) 



T = 



where the generalized quadrature = (ae - *^ — 

a jf e lc ^)/2i. Then, position and momentum operators are 
measured when choosing (f> = —n/2 and = 0, respec- 
tively. To apply this technique, it is required a particular 
initial state of the internal states, so the measurement of 
(x(t)} would have to be done in two parts, one associated 
with each of the suitably projected states, |+^) and |— 

We turn now to the possible simulation of Klein's para- 
dox. In 1929, Klein noticed [l|| the anomalous behavior 
of Dirac particles in regions where a high potential energy 
V exists: = h£ ) +V1 2 - When V > 2mc 2 , negative- 
energy electrons (components) may swallow V, acquiring 
positive energy and behaving as ordinary electrons, while 
leaving a hole in the Dirac sea. This stems from the fact 
that the relativistic energy related to Hy may be recast 
into (pc) 2 = (Ey^ -V+ mc 2 ){E^ ] - V - mc 2 ), which is 
positive when either both factors are positive or negative. 



In the second case, the total energy E^p < 



-mc 



V 



can be larger than mc , as noticed by Klein. In this 
case an e~ — e + (electron-positron) pair could be cre- 
ated from V. The sudden raise of the constant poten- 
tial V(|a)(a| + at a certain time t = to after an 
evolution associated with could simulate this phe- 
nomenon. We suggest to produce potential V with the 
required characteristics through a fast and homogeneous 
Stark shift in both internal levels. The natural way to de- 
tect Klein's paradox, assuming an initial positive-energy 
internal state \a) (po = 0), is via measurement of nonzero 
population in the negative-energy component \b). 

Finally, in as much as massless chiral fermions in 
condensed matter, we could also produce a 1+1 axial 
anomaly [ItJ by changing Eq. ([TO]) into Hp = ca x p x + 
qEx, where x cx (a + a 1 ) is a motional displacement. 

The basic ingredients to implement the ID Dirac dy- 
namics of Eq. (|10| with a single trapped ion are two in- 
dependent electronic (internal) states coupled to a one- 
dimensional motional degree of freedom. The required 
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FIG. 1: Relevant energy levels (not to scale) of a 25 Mg + ion. 
Shown are the ground-state hyperfine levels supplying the two 
internal states, |o) and \b), and the equidistant harmonic os- 
cillator levels related to the harmonic axial confinement in 
a trap similar to that described in [lH [l9l |. We subsumed 
excited levels of the P\/2 and P3/2 states. We depict the 
resonant transition state sensitive detection and the relevant 
types of two photon stimulated Raman transitions between 
states \a,n — 1} and \b,n — 0,1,2): red sideband (dotted 
line), carrier (solid line), and blue sideband (dashed line). 



states could be composed by two ground-state hyperfine 
levels of an earth alkaline atomic ion, e.g. of 25 Mg + by 
|F = 3; m/ = -3) and \F = 2;m f = -2), |o> and |6), 
respectively, as depicted in Fig. 1. A constant external 
magnetic field will define the quantization axis and lift 
the degeneracy of levels being potentially useful to pro- 
vide additional states (like |c) and \d) in Fig. 1) necessary 
for higher dimensional simulations, see Eq. @. 

At the start, the ions will be laser cooled close to the 
motional ground state and optically pumped into state 
I a) [2(|. One red/blue sideband and one carrier tran- 
sition will be driven simultaneously (lH 21 1 to imple- 



ment the desired dynamics of Eq. ([5]) via two-photon 
stimulated Raman transitions. To measure the ion po- 
sition we rely on the mapping of motional information 
on the internal degrees of freedom 14, [l5| and take 



advantage of the high fidelity of state sensitive detec- 
tion realized by an additional laser beam, tuned to a 
cycling transition [23], coupling state \a) resonantly to 
the P3/2 level. Considering the available laser intensi- 



ties [23(1 , all necessary Raman beams could be derived 
from a single laser source. We split the original beam 
and provide the necessary frequency offsets, phase control 
and switching via multi-passing through Acusto-Optical- 



Modulators [2l|.The number of laser beams could be fur- 
ther reduced by electro optical modulators to provide red 
and blue sidebands simultaneously [TJ [l2j]. To imple- 
ment an overall shift in the potential (Stark shift of level 
I a) and \b) in Klein's paradox) without changing their 
mutual energy difference we chose the directions and po- 
larizations of the Raman beams appropriately. 

In conclusion, we have shown how to simulate the Dirac 
equation in 3+1 dimensions for a free spin-1/2 particle, 
and quantum-relativistic effects like the Zitterbewegung 
and Klein's paradox, in a single trapped ion. We have 
studied the 1+1 case, where experimental needs are min- 
imal while keeping most striking predictions. We be- 
lieve that these simulations open attractive avenues and 
a fruitful dialog between different scientific communities. 
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